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^ ; Abstract 

^ . ' The derivation of the angular spectrum of temperature perturbations of 

I ' the cosmic microwave background reUes on the quantization of field and metric 

^J^i perturbations in the inflationary phase. The quantization procedure thus de- 
serves a close examination. As the background spacetime on which the degrees 

' of freedom that are quantized live is curved, the methods of quantum field 

^> . theory in curved spacetimes are applicable. Furthermore the dynamic system 

I that is quantized contains a constraint, which adds an interesting problem to 

' the quantization of the system. This article investigates the treatment of the 

■ constraints and the Weyl quantization of the system. A preferred quantiza- 

CSl ! tion procedure is identified, which renders the field and metric perturbations 

' non-local. 
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1 Introduction 

The theoretic framework of hnear cosmological perturbation theory has been used to 
predict the angular spectrum of temperature perturbations of the cosmic microwave 
background with great success. As the derivation of the perturbation spectrum rehes 
on the assumption of an inflationary phase in the early universe, the success of the 
predictions provides a cornerstone argument for the assumption of an inflationary 
phase in the early universe. 

The derivation of the perturbation spectrum starts out with a generic scalar field 
setting for inflation where a semiclassical treatment is employed. The semiclassical 
Einstein equation, see [21] for an early discussion, provides an approach which is 
widely used in the context of quantum field theory on curved spacetime. In this con- 
text, matter is assumed to be modeled by quantum field theory, while the spacetime 
geometry is treated classically. The matter contribution is then represented by the 
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expectation value of the energy-momentum operator. In the case at hand a different 
setting is chosen, where both the metric and the field are split up into a classical 
background part and a small quantum perturbation, as e.g. in [1], [TT] and [27] . 

As a rationale for this kind of semiclassical approach one can point to the very high 
energy density that is expected to govern the inflationary phase suggesting the need 
to account for quantum effects of gravity. However, without a generally accepted 
theory of quantum gravity any quantization scheme for the metric perturbation 
has to be considered ad hoc. The question concerning validity of the quantization 
procedure based on the perturbative splitting will not be pursued in the present 
work. 

The quantum state of the perturbations during inflation is usually assumed to be 
close to the Bunch-Davies state [6j on de Sitter spacetime. The infrared asymptotics 
of the two point distribution is then used in the derivation of the angular spectrum of 
perturbations. The choice of state can be justified by invoking the infrared stability 
of quantum field theories on de Sitter spacetime. The validity of this concept for a 
general set of models is another interesting topic, that will not be explored in the 
present article. For recent results and an abundance of references, see PTS] and |25j . 
There have also been several investigations of the influence of the initial state on 
the perturbation spectrum which found constraints on the sets of possible initial 
states from the measured perturbation spectrum. For an overview of the different 
approaches that have been pursued, see [5], [20], [22] and references therein. 

It should be noted that in the course of the derivation, the quantum two point 
distribution is tacitly identified as a classical correlation. This feature was consid- 
ered worrysome by |26j who proposed an alternative derivation of the spectrum of 
perturbations, which included an explicit transition from the quantum state to an 
effectively classical situation. Although the proposal made by [^H] seems somewhat 
sketchy and incomplete, we contend that the quantum-classical transition deserves 
some further investigation. The present article will however not be concerned with 
this question. 

The derivation of the angular spectrum of temperature fluctuations from the in- 
frared asymptotics of the two point distribution interpreted as a classical correlation 
is purely classical. As we are interested in the quantum aspects of the problem, we 
will not investigate this part of the derivation. 

In the present article the quantization of the dynamic system of perturbations 
during inflation will be investigated. As preparation, section [2] introduces the model 
that is investigated in the following. We then go on to prove a no-go result concern- 
ing simultaneous canonical quantization of field and metric perturbations in section 
131 Section m is dedicated to the investigation of the dynamic system of perturbations 
in a phase space setting using standard tools for dynamic systems with constraints. 
The Dirac bracket governing the classical dynamics is identified and in section \5\ 
an argument is proposed, which enables a unique choice of dynamic variables. The 
preferred quantization procedure implies non-locality of the field and metric pertur- 
bations which provides an interesting connection to non-commutative inflation [T]. 
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Lastly, Weyl quantization of the dynamic system with constraints is discussed in 
section O 

2 Linear Cosmological Perturbation Theory 

Cosmological perturbation theory is a special case of perturbation theory in general 
relativity. The present article will only be concerned with linear perturbation theory, 
which means that only the tree level and the first order in the perturbations are con- 
sidered and higher orders are neglected. The basic idea is to split all quantities up 
into a background part and a perturbation. In the case of geometric quantities, most 
prominently the metric, this is nontrivial due to the role of diffeomorphism covari- 
ance in general relativity. It turns out that only part of this covariance survives in 
the setting of linear perturbation theory, which is interpreted as covariance under in- 
finitesimal diffeomorphisms and labeled gauge freedom. For physically unambiguous 
results it is reasonable to define gauge invariant quantities. 

Using a 3 + 1 decomposition of the background metric, as reviewed in [30j, it is 
possible to classify perturbations by their transformation behavior under symme- 
tries of the spatial section. This leads to the definition of scalar, vector and tensor 
perturbations. In the present article we will only be concerned with scalar pertur- 
bations, thus vector and tensor perturbations will be dropped, after they have been 
identified. 

The treatment in this section mostly follows [27] where no other sources are cited. 

The background spacetime is assumed to be have a Friedmann- Robertson- Walker 
metric with flat spatial sections with line element 

ds^ = a^{T]) {dif — dx^^ . 

Here, rj is the conformal time and the spacetime is obviously conformally fiat. For 
some given quantity X the derivative with respect to conformal time is denoted as 
X' . We also denote % = ^. The time dependence of the scale factor a depends 
crucially on the matter model coupled to the spacetime via the Einstein equations. 

As preparation for the formulation of the considered model, a short introduction 
to linear perturbation theory is in order. The approach taken here differs slightly 
from the usual treatment in the literature. The presentation here is slightly differ- 
ent in style from that of |1] and ^Tj, but clearly leads to the same results. We 
thank A. Reichold for pointing out this quite appealing formulation of the founda- 
tions of linear perturbation theory and gauge transformations. The key element for 
our description of linear perturbations and gauge freedom is the interpretation of 
perturbed quantities as one-parameter families of tensor fields. 
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Definition 2.1. 

• Let O C M an open neighborhood of and let T : (9 — ^ r(TJA^) a smooth 
one-parameter family of smooth r-s-tensor fields. Then Tq := T(0) is called the 
background quantity and 5T := -^T[\)\^^q\ is called the linear perturbation. 

• Let $ : (9 — ^ DiffTW a smooth one-parameter family of diffeomorphisms gen- 
erated by the vector field X, i.e. $(0) = id>i, X = ^^'(A)|;^^q and for the 

pull-back (T;$)(0) = idr(TjAi)- Then we define df := |^((T;$)(A)T(A)) l^^gA 
as the transformed perturbation while the background quantity is unchanged, 

fo-.^ (r;$)(o)To = ro. 

The above definition has the benefit that it is mathematically concise and yields 
the transformation law of the perturbations under gauge transformations in a very 
transparent way, 

5f = 5T+ X^xTo =: 5T + ^(Tq. 

In this equation ^ = XX is often interpreted as an "infinitesimal vector field" to 
justify the linearization. In any case the linearization is seen to be essentially a 
linearization of the Taylor expansion in A of one-parameter families of tensor fields. 
Using the Levi-Civita connection and local coordinates with respect to the back- 
ground metric go we get 

n=l n=l 

where the multiindices ct„ and p„ are equal to the multiindices cr and p respectively, 
except that the index at the n-th position, i.e. an or pn respectively, is replaced 
by a. Due to the metric compatibility of the Levi-Civita connection this yields an 
especially simple transformation law for the metric perturbation 5g given by 

It is customary to classify metric perturbations using a 3 -|- 1 decomposition of 
the background spacetime into spatial section and time direction if the background 
metric is sufficiently symmetric. In our case the background metric is a Robertson- 
Walker metric with fiat spatial sections, therefore such a decomposition is indeed 
beneficial. The perturbations are decomposed into scalar, vector and tensor part, 
dependent on their transformation behavior under the symmetry group of the spatial 
section. 

Let 7 denote the metric of the spatial section E and let X\i denote the application 
of the Levi-Civita connection of (E, 7) on X expressed in the coordinates of the 3-1-1 

decomposition. In this context, Sqqq =: Si is considered scalar; Sgoi =: S2\i + Vn is 
decomposed into a pure gradient scalar part S2 and a divergence free vector part 
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Vi with vi\i — 0; Sgij —: Ss\ij + s^^ij + V2i\j + V2j\i + tij is decomposed into two 
scalar parts, one pure gradient S3 and the trace S4, the gradient of a divergence free 
vector part V2 with ^2*11 = and a trace and divergence free tensor part Uj with 

t i t j^i 0. 

Counting the degrees of freedom, one naively gets 4 scalar, 4 vector and 2 tensor 
degrees of freedom. However, the vector fields { (generating gauge transformations) 
can be decomposed into a scalar part ^5 = ^°9o + and a vector part — x^di 
with = 0, with two components each. This means that the gauge freedom reduces 
the degrees of freedom such that 2 gauge invariant degrees of freedom of every type 
remain. An important point of this decomposition is that the different types of 
perturbations do not mix under gauge transformations. 

In the following the necessary formulae for a basic treatment of cosmological ques- 
tions in the context of linear perturbation theory are derived. As explained above, 
the background metric used for cosmological perturbation theory is a Robertson- 
Walker metric with fiat spatial section. We restrict now to scalar perturbations, so 
the perturbed metric can be expressed using conformal time as 

ds'^ = a^{ri) [{I + 2A)drf + 2B\idr]dx' - ((1 + 2D)5ij + 2E\i^)dx'dx^ 

where A, B, D and E are functions of all coordinates. It should be noted that 
A, B, D and E are not exactly the same as the s„ used above, but differ by a 
factor of 2o?{rj). This yields slightly different behavior under gauge transformations. 
As we restrict to scalar perturbations we also restrict to scalar gauge vector fields 
^ — ^°9o -I- This yields transformation laws 

A=A + + 5 =5 -t- C° - 



Obviously any two of the perturbations can be brought to vanish simultaneously ex- 
cept A and D, as they are both only changed due to This means if A and D can 
be brought to vanish simultaneously in some gauge, this has a physical implication. 

An often used gauge is the conformal gauge, where B^ = E^ = 0. To separate phys- 
ical contents from gauge effects, it is useful to work with gauge invariant variables. 
As we have four scalar perturbations and two gauge parameters we can derive two 
gauge invariant scalar perturbations. The ones usually defined in the literature are 
the Bardeen potentials 

* =^ - ^ ((S - E'))' <^>^D-n{B- E') 



As matter model, we use a simple scalar field with a so far unspecified self- 
interaction potential V . The field satisfies the Klein-Gordon field equation 

dV 

□0-- — (0) 
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where V{(f)) is some potential. The energy-momentum tensor takes the form 



In keeping with the idea of hnear perturbation theory also the scalar field is 
decomposed into a background part and a perturbation = 0o + <^0- The gauge 
transformation behavior of the field perturbation is 



Again in conformal gauge 6(f)c = X'l the following the subscript indicating confor- 
mal gauge will be suppressed and whenever expressions contain 50 they are to be 
understood in conformal gauge. 

To formulate the Einstein equation it is necessary to calculate the Einstein tensor, 
which is not gauge invariant. However G^— SvrGTjJf can be cast into a gauge invariant 
form, as it vanishes in any gauge. This implies that the Einstein equations can be 
formulated using only gauge invariant variables. This means we can calculate the 
Einstein tensor and the energy momentum tensor in some arbitrary gauge, and the 
resulting Einstein equations will be true for any gauge. Hence, we can replace the 
perturbation variables in the Einstein equations by gauge invariant variables that 
are equal to them in the chosen gauge. The resulting equation, containing only 
gauge invariant variables, still holds in any gauge. A calculation using the above 
strategy is conducted in [27] and we will only give the results here. 

The idea of perturbation theory is to assume the equations to hold order by order. 
This means we will have a background Einstein and Klein- Gordon equations and 
perturbations to these equations which will be assumed to hold separately. The 
background Einstein equations are 



The perturbation Klein-Gordon and Einstein equations will lead to a set of four 
independent equations. From the Einstein equations we get 



50 = 50 + eVo 
which implies the following gauge invariant perturbation 

X = 5ct>-<P',{B + E') 



and the background Klein-Gordon equation reduces to 

0'o' + 2H0'o + K*(0o)a' = 



(2.1) 



(2.2) 



^" - + 2 f?^ - ^ 1 ^' + 2 ( - 1 ^ = 



(2.3) 



(2.4) 
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and the first order of tlie Klein-Gordon equations yields 

X" - Ax + 2W + V,^^ (0o)«'x = 40[,*' - 2V,^ (0o)«'^ (2.5) 



Tlie coupled differential equations (12 ■4p . fl2.3l] and fl2.5p govern the dynamics of 
the gauge invariant perturbations ^ and x- The present article will investigate such 
perturbations using these equations and $ will not be treated as an independent 
variable due to equation (12. 2p . An alternative derivation of the dynamic system 
is possible by means of a Lagrangian approach reviewed in The Lagrangian 
approach will be investigated in section HI 



3 A no-go Result 

In this section it will be shown that \i/ and x cannot both satisfy canonical commu- 
tation relations. To show this, the dynamical system consisting of equations (12. 3p 
and (12. 4p is used. The first objective is to simplify and unify the notation of the 
equations. 

For simplified and homogeneous display of the equations, the following calcula- 
tions use the parameter z = ^ and the slow-roll parameters 



. = 1+.-^ (3.2) 

With the additional use of of the Klein-Gordon equation (12. ip for the background 
field the following useful relations can be derived. 

-=(l + e-5)H — = (2 + 26-36 + 26^ - 3e6 + 6^-^) 

^=(1-^)^ ^=C2-^e + 2eb)n^ 

#=(1-^)^ ^ = (2 -£-3<5 + £(5 + (52- 1^^7/2 

00 00 V ^/ 

V,^ = - (3 - V,^^ a" = ({3 - S){e + S) + Q 



- =2{e - 5)n 
e 
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These relations can be used to unify the notation for the dynamic system such 
that only a, Ti, 6, e and z appear in the coefficients. 



-X 



Thus we can derive 



X 



^ T/ ^ T 



eT-ia 



ea 



(3.3) 
(3.4) 

(3.5) 



from equation (13. 4j) . Taking the time derivative of equation (13. 4p and applying 
equation 03.31) one gets 



+ (1 - 25)H^' + {l + e- 26)n^^ 



ix' + (1 - 6)nx) (3.6) 



Now, application of equation (13.51) to eliminate x yields 

x' = --^' + ^(A + {e-s)n')^ 



(3.7) 



After these preparations we can now prove the announced no-go result. We em- 
ploy symbolic notation in terms of integral kernels of distributions to avoid the 
technicalities of defining "equal time" in the distributional sense. 

Lemma 3.1. 

If there is a homogeneous bidistribution with integral kernel f{x—y) and the following 
relations hold 



[^(x),^(y)]| =[vl/'(x),vl/'(j/)] 



x"=y" 



-- 

f{x-y) 



then 



[xix),xiy)] 



[xix),x'iy)] 



[x'{x),x'{y)] 

z^^A^' + en 



[^{x),^'iy)] 



(3.8) 
(3.9) 

(3.10) 
(3.11) 



Note that a, e, H and z only depend on time and the expressions investigated here 
are evaluated at equal times. Therefore their time dependence is suppressed. 
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Proof. 

We use equations fl3.5p and 03.71] to compute the commutation relations 



[x'{x),x'iy)] 



'(A" + {e- 5)%'') [/\y + (e - 5)7/' 



{ea]m 



[x{x),x!{y)] 



(eHa) 



Sz"" 



(ea) 



[^{x),^'{y)] 



(May ^ '^W-*' 



This proves the claim. 



□ 



As a direct corollary it is not possible for \1' and x to both satisfy canonical 
commutation relations (CCR). Likewise they cannot be in a Hadamard state at the 
same time. However, these states are assumed to contain all physically acceptable 
states for a quantum field theory on a curved spacetime. An early work pointing out 
the importance of Hadamard states is for a recent application see e.g. [8j. The 
no-go statement is a consequence of the fact that the constraint equation reduces 
the number of independent degrees of freedom to one, so the result is not surprising 
per se. It is interesting to derive a dynamic equation for x alone from 02.51) using 
equation 03.51) and 



A$ + £7{2^ = ^(x' + 5-Hx) 



(3.12) 
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which can be derived from equation (13.61) . The resultant equation 



;A + eH^)-(" + 2H{A + e5H^)-(' + 



-A^ + (-2e + 3(5 - £5 - 52 + 1 J 




e + 5-e5 + 5'^ + — \eH^\x = Q (3.13) 



is not a normal Klein-Gordon equation but a dynamic equation of fourth order. 
Thus it is not to be expected that the CCR or the Hadamard condition are of any 
significance for x- The field ^ however fulfills a Klein-Gordon type equation, so at 
first glance it seems reasonable to assume standard CCR for this field as a means 
of quantization and treating x as a derived field whose commutation relations as 
derived from the CCR for \l/ bear no significance. However, this naive approach will 
be proven misguided by the following investigation. 



4 Canonical Symplectic Form 

The usual choice of the field to quantize is m = z'^ + ax which is justified by the fact 
that the quantity IZ = —'^ can be interpreted as gauge invariant curvature perturba- 
tion, which is the quantity that is used to derive the temperature perturbation. In 
the present treatment we take the point of view that the physical significance of the 
quantity does not justify choosing it as canonical variable for quantization. There- 
fore we are looking for a structural argument to identify a preferred quantization 
procedure. 

In this section the dynamic system of the spatial Fourier modes is reviewed from 
the angle of symplectic geometry. To this avail, the dynamical system has to be cast 
in Hamiltonian form, which implies the identification of the Darboux coordinates 
of the symplectic form. To achieve this, it is necessary to use the action presented 
e.g. in and to derive the corresponding Hamiltonian system. Then the impact 
of the constraint on the symplectic structure is discussed. The following discussion 
will employ Hamiltonian and Lagrangian densities, where the term "density" will 
always be suppressed. 

For the treatment of a dynamical system with second class constraints several 
methods exist. The dynamical variables, including the Lagrange multipliers, may 
be combined into a reduced set of dynamical variables in such a way that the result- 
ing Lagrangian contains no constraints and only depends on the new variables. This 
is the strategy which is pursued in the standard literature on cosmological pertur- 
bations during inflation, especially [1]. However, finding an appropriate reduced set 
of dynamical variables is not a straightforward task and is usually done by simply 
choosing the appropriate ansatz to begin with. As the resultant dynamical system 
is a generic Hamiltonian system, quantization is straightforward. 

A second method consists in the straightforward construction of a Hamiltonian 
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containing the Lagrange multipliers and the subsequent application of a standard 
procedure to eliminate the Lagrangian multipliers. This procedure does not lead to 
an obvious reduction of phase space, however it can be seen that the phase space 
is not symplectic and the complement of the maximally symplectic subspace has 
a dimension which is identical to the number of constraints. This procedure has 
the disadvantage that the geometry of the phase space is unclear and therefore 
quantization of the dynamical system is not straightforward. 

Yet another method consists in ignoring the constraints at first, to produce an "un- 
constrained Hamiltonian" with a corresponding canonical symplectic form. Then the 
constraints are used to derive the so-called Dirac bracket from the Poisson bracket 
of the unconstrained system. The tensor corresponding to the Dirac bracket will be 
called the "Dirac form" in the present treatment. The Dirac form is degenerate in 
a number of dimensions corresponding to the number of constraints and symplectic 
in the remaining dimensions. The present treatment will pursue the derivation of 
the Dirac bracket to shed some light on the quantization of the system. The classic 
reference for this method is [TU] but it is treated to varying degree of detail in several 
textbooks, e.g. 0, [IH| and [S2|. 

We would like to point to [7J, where a similar scenario is investigated but vector 
and tensor perturbations are also taken into account. However, that work applies a 
gauge fixing procedure to achieve a reduced phase space, which leads to a different 
treatment than pursued in the present work. 

As a first step it is convenient to rescale the fields \& and x so they have the same 
dimension. Therefore we define the rescaled fields 



The Lagrangian given in [1], equation (10.68) is given here in terms of the rescaled 
variables, suppressing the total divergences and simplifying by setting $ = — \E' and 
L := 4 A (5 — E') (note the different notation). From the treatment in [4] it is clear 
that the terms introduced by switching to gauge invariant variables cancel out. 



The fields are treated as real fields here, but a generalization to complex fields is 
st raight forward . 

The next step in the treatment is to consider the "unconstrained Lagrangian", 
thus setting L = 0. We will additionally switch to a "mode Lagrangian", using 



i/j := z'^ 



X := ax 



1 r 



3(V'')^ + 6(e - 5)n'ip'il^ +{e- 3e^ + 6e5 
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modes tpi, := J e^^^tj^^x). The modes only depend on k = \k\ due to the fact that 
the dynamics are spatially homogeneous and isotropic. 



2e 



-3(V'L)' + 6(5 - 5) WfcVfc + [(^ - + Qe5 - - k^]il^l 



1 

+ 2 



(A;)2 - 2-HA;Afc + 



1 - 3£ - 35 + e5 + 52 _ ^ j -H^ - A;2 



+ 4?^^/'lAfc + 2(1 - 25 + 5)n''xl)kXk 

Note that the mode Lagrangian is not the Fourier transform of the Lagrangian, 
as a Fourier transformation would not map products of functions to products of 
their Fourier transforms but to convolutions. The spatial integral of the Lagrangian 
coincides with the mode integral of the mode Lagrangian, which means they lead to 
the same Lagrangian function and thus to the same dynamics^ Due to the fact that 
the dynamic equations do not mix modes, no mode-mixing occurs in the Lagrangian. 

From the mode Lagrangian one calculates the canonical momentum modes 



-V. + 3— 



6 



-gHfc + {e 



<^ Al, — Kfc + HXk 



which allows the derivation of the unconstrained Hamiltonian for the modes 

1 



0,fc 



2e 



+ ^^fc + "H/tfeAfe + 



-7e + 96- 3e6 - 36^ 



3^) + 



+ ^enUkXk-2{i-s)n^i^,Xk 



The unconstrained phase space of Fourier transformed initial values is assumed 
to be (^(M+)) in the following, which means that initial values are Schwartz 



functions in the spatial coordinates. The Hamiltonian density at fixed time is then 
a smooth functional : (^y{R+)Y — > y{R+). The initial value phase space 



^This is essentially an application of the Fourier-Plancherel theorem. 
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for a single mode is because the modes are simple functions of time. The mode 
Hamiltonian is thus a smooth function Mo^k '■ — ^ l^- 

Any smooth phase space "density" ^ : (^^(M"*")) — > ^(R+) can be integrated 

to yield a smooth phase space "function" F : (^^(M+)^ — > M. Furthermore, for 

a polynomial / : M+ x M with V/c G M+ : /(A;, 0) = 0, a smooth density ^ 

can be constructed as 

(91, 92, 93, 94) (k) ■■= fk(9i{k),92{k),93{k),94{k)) (4.1) 

where ^[91,92,93,94) e <5^(ffi+) if 9i, 92, 93, 94 e <5^(IR+). Let us call the space 

of such polynomial densities r'i>oi,4(^J^(M+)^ and the functions derived from such 

densities by integration as J^Poi,4(-5^(K+)V As such polynomials encompass the 
mode Hamiltonian and the field modes, this suffices for our purpose. It should be 
noted that for a general smooth function / : xR'' ^ M with VA; e R+ : f{k, 0) = 

a smooth function F : (^y(M.'^)) — > R not necessarily exists. 
Now consider the Poisson bracket on a single mode 

{, }k : C°°(m^ ^ m) X C"^(m^ ^ m) — >R (4.2) 

The Poisson bracket for a phase space of several modes is always diagonal in the 
modes due to the symmetry of the dynamic system. It is obvious that this bracket 
can be lifted to densities created from polynomials as 

{, }v : Vp,,,,{y{R+)) X Ppoi,4(^(R+)) =i^(K+) 

{^,^}j,{k):^{nk)Mk)}k (4.3) 
and thus a Poisson bracket on J-poi4(<^(R^) ) is defined as 



{F, G} = J '^}v{k)47redk = J {^k, %}k4Mk (4.4) 


Note that we use the three dimensional volume element Ank'^dk to account for the 
fact that k is actually a three dimensional vector. Only one integral is performed 
because the mode Poisson bracket is diagonal in the modes. 

As the canonical variables of the unconstrained Hamiltonian are the Darboux 
coordinates of the canonical symplectic form on the unconstrained phase space, the 
symplectic form corresponding to the mode Poisson bracket takes the standard form 



O2 I2 



-I2 O2 J 

on the linear space of first derivatives of mode functions '■ 
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The constraint equation imphes a projection to a subspace of the phase space. 
Contrary to the projection to an energy hypersurface this projection constrains the 
accessible phase space for all solutions of the dynamic system so it can be interpreted 
as a genuine reduction of phase space. In canonical variables the primary constraint 
derived from equation 03. 5p reads 

Ci = Hfe - H\k = (4.5) 

and a secondary constraint can be derived as 

C2 ■■= -^{Ci, ^0. J - = f^k- + Sm, = (4.6) 

using the first constraint. The tertiary constraint vanishes 

dC2 

Cs = {C2, ^o,k} + = 

applying the first two constraints, so there are exactly two geometrically independent 
constraints, which means that the phase space dimension is reduced by two. It is 
worth emphasizing that these constraints are satisfied by all valid solutions of the 
dynamic system, so all solutions of the mode dynamic system can be described as 
living on a two dimensional reduced phase space. 

For the treatment of constraints it is of great importance whether they are first 
or second class constraints. In the case discussed here, only second class constraints 
occur, due to the following lemma. 

Lemma 4.1. 

The constraints fl4.5p and f l4.6p are second class constraints for k 0. For k ^ 
the Dirac form diverges no stronger than of order 

Proof. 

To fix notation let {Mab)a,b denote the matrix with components Mat- If the con- 
straint matrix (Cab)ab '■= i{Ca,Cb}) is invertible, the constraints are second class 
constraints. Using the definition of the Poisson bracket 

df dg df dg df dg df dg 

if, 9} ■= TT^TFT + 



dip^dllk dXkdnk dUkdipk dukdXk 
one gets 

Cu = {CuC2} = -n = 

which yields, due to antisymmetry 

{Cab)a,b = ^ ( 1 ) 

which is invertible as long as 7^ 0. As the inverse diverges as also the second 
part of the claim is true. □ 
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First class constraints usually occur in the context of gauge theories so one might 
expect that a full treatment of scalar perturbation theory, not using gauge invari- 
ant variables, might be significantly more complicated, especially when it comes to 
quantization. Such a more extensive investigation is not attempted in the present 
work. 

In the present case, the constraints are first class constraints for the zero mode 
k = 0. This is however an artifact of the mode setting and it will be seen that a 
complete formulation of the Dirac bracket on the full phase space of all modes is 
nevertheless possible. 



The Dirac bracket is derived from the Poisson bracket and the constraint matrix 



as 



{f,9}D■■={f,9}-{f,C^}iC'%{C„g} 
and a straightforward calculation yields the mode Dirac form 



/ 



en 






V 



en 

e2f 

esnl 

fc2 



en^ 






esn^ 



1 - 



fc2 

en? 

fc2 

at en^ 







(4.7) 



which obviously diverges as k""^ for k 0. The Dirac bracket on all modes is 
defined analogously to the Poisson bracket on all modes. On the linear space of first 
derivatives, we cast it in the symbolic form 



V h ) 



(91 \ 

92 
93 
\94 / 



i( flik) \ 
f2{k) 

h{k) 

V U{k) J 



( 9i{k) \ 
92{k) 
9zik) 
\ 940) J J 



d^k 



(4. 



D,k 



where all the functions fi and Qi, symbolizing the values of first derivatives of smooth 
densities ^ : (^^{M.'^)^ — > are Schwartz functions in k. As the mode 

Dirac form given in equation (14. 7|) behaves like for A; — j- and the volume element 
of the integral in equation (14. 8 p gives a factor of k"^ the limit of the integrand for 
— is finite. Thus the integrand is a Schwartz function and the integral is well 
defined. 

As a last step, the constraints can be factored out of the phase space, such that 
the resultant factor space can serve as physical phase space. This procedure will be 
investigated to some more depth in the context of quantization. 



5 Preferred Dynamic Variables 

We now turn to investigating which canonical variables are suited for quantization 
of the dynamic system with constraint. The canonical variables can be chosen freely 
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and only need to satisfy the requirement {X,P}d = 1. This freedom imphes the 
apphcation of the constraints as well as symplectomorphisms. To make a connection 
to sectional the possible canonically conjugate variables P for X = Uk = if^k + -^fc 
and for X = \il^k be identified here. 

Considering X = the corresponding conjugate variable has to satisfy 
eH dP eU dP sH'^ dP ^ - sH'^ - e6H^ dP 

+ ^^7^ + 



/c^ dt/jj^ k"^ d\k k"^ dllk k"^ 
which is indeed fulfilled by 

4 



{e - 5)ni^k + (i + 3^) - |nfe + Kk 



which implies that the standard quantization procedure is indeed justified. 
For X = the corresponding conjugate variable has to satisfy 

eHdP_ eW^ dP -e6n^ dP _ ^ 
zk"^ dXk zk"^ dUk zk'^ dnk 

which can be fulfilled by 

P = aX' = V'fc + 7T Afc - —Uk 

z 3 z 3z 

namely if a = This result is peculiar in that it includes a prefactor of fc^. 

A comparison with equation (15.181) shows that the treatment of this subsection 
reproduces the result of the previous subsection with respect to the relation of the 
commutators (or in the classical setting Dirac brackets). However, the treatment in 
the present subsection allows to take the point of view that the standard treatment of 
quantizing u is canonical in the sense of Dirac quantization, while the straightforward 
CCR quantization of \1/ is not. 

The choice of the field u as canonical variable can be justified by rather natural 
technical requirements, as formulated in the following result. 

Theorem 5.1. 

Let A and B be smooth real-valued functions of time hut independent of k and let 
B he strictly positive. Assume X^ = Ai/)j^ + BX^ then B = 1 and A = 1 or A = 
1 + '^'"^ ihe only solutions compatible the technical requirements 

(a) Pk = X', with{Xk,Pk}D = l 

(b) the dynamic equation for X is a differential equation of at most second order. 
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Proof. 

The first requirement leads to the equation 



-AB + 



-A' 



(5.1) 



^ eH{AB' - A'B) + {k'^ -eH^)B'^ + {l-e)eH^AB + e^H^A^ = k'^. (5.2) 

Because B does not depend on k, B'^ — 1 and thus B — 1. This yields 

A' - sHA^ - (1 - e)nA + A' - (A-l)(Ae + l)n = (5.3) 

which can be cast in a simpler form using D := A — 1 

D' - D(De + (1 + s))n ^ (5.4) 

such that one obtains a Bernoulli differential equation. This type of equation is 
solvable with trivial solution D — and a nontrivial solution which can be found by 
using E — ^ leading to 



E' + {1 + e)nE ^ -en 
a readily solvable linear equation. The solution for A is 



1 



n 



J e{ri')a^{ri') drj' + K 



(5.5) 



(5.6) 



Wo 



where K and r]Q account for one free constant. 
Now abbreviating p :— A — 1 one readily sees 

p' = eHp'^ + (1 + e)Hp 

p" = 2e^n^p^ + 2e{l + e- 5)V?p' + 2(1 + £ + _ s5)V?p 



(5.7) 
(5.8) 



Defining Up^k = Uk,o + pV'fc ^hc dynamic equation can be derived from 

u';, = u'^ + p"il, + 2^'^' + pV'" (5.9) 
using the dynamic equations 



2 + 2£ - 35 + 2£^ - 3£(5 + - ^ ) - 



xl^l = {2 + 2e- 46)71^', + 
and the constraint 



6'^ 



2s - 5 + 2e^ - eS - 5^ - ^]n'' - e 

n J 



(5.10) 



V'. (5.11) 



(5.12) 
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It takes the form 



u' 



2e'p' + 4(1 + e - 5)ep + 2 + 2e - ^6 + 2e' - ^e5 + 5'' - ^ \ - k"^ 



p,k 



Uph 



+ 2 



- e{l + e)p+{e-5){l + e- 5) n^pt/^k 



(5.13) 



As the secondary constraint can be cast in the form 

^k = -^{- «M + (^P + 1 + ^ - mup,k) (5.14) 

one readily sees that the coefficient of i/?^, in equation (15.131) must vanish to satisfy 
the third requirement, as orders of k translate to spatial derivatives. The condition 
—e{l + e)p'^ + [e — 6){1 + e — 5)p = implies the claim. □ 

Although M = i/j + A is not uniquely determined by theorem 15.11 it is clear that 
the alternative solution implies restrictions on e and S. In this sense, u is the unique 
suitable choice for all background configurations. We will point out one feature of the 
alternative solution here which is especially unwanted in the context of inflationary 
models. 



Lemma 5.2. 

-h.^ . - (£: 



If 6 = 0{e), then p = implies {e - 6)' = {e - 6)H + 0{e^). 



Proof. 

Using relation fl5.7p one can derive 



6 1+6 e — 5 1+e — S p l + e 



which can be recast using 7 := e — 5 to yield 



(5.15) 



_ 7(1 + 7)(7(1 + 7) + (1 + ^)^ + 27(1 + 2.)) 
^ " (l + e)(l + 27) ^ ^^-^^^ 

The claim follows by Taylor expansion of this expression. □ 



If one invokes the slow roll assumption and assumes the higher order terms to be 
small, one is left with a readily solvable differential equation, that yields e — 7 oc a. 
Thus, contrary to the assumption, e and 7 are not small and slowly varying, which 
shows that the restrictions on e and 6 imposed by the alternative solution are incom- 
patible with slow roll inflationary models. We will therefore discard this solution 
and regard X = 1/? + A as a preferred solution. This result gives a justification, from 
a more structural point of view, for the usual choice of canonical variable in the 
standard literature on cosmology. 
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To generalize the result presented here, one might relax the first requirement in 
theorem I5.H by allowing P^. = DXj, where D is some smooth function of time, but 
does not depend on k. One would expect that this introduces at least a freedom of 
the form X = + A) for the preferred solution where F is some smooth function 
of time. 

The standard approach to quantization assuming OCR for the field u = X + ip = 
ax + z"^ was shown to be justified. It is thus interesting to investigate what the 
CCR quantization of u implies for the commutation relations of \E' and x- 

Lemma 5.3. 

// standard CCR for u are imposed, ^ and x o'^e non-local. 



Proof. 

As a first step we assume 



0. 



Using equations 03.51) and 03.71) the following relations can be derived 

u =— ^ H z-^ 

en e 

u' =^/\^ + {l + e-5)Hu 
BrL 



(5.17) 



[u{x),u'{y)] 



z 



^0 = yO ^2^2 

Now the last relation can be recast in the form 

x"=y' 



? z 

^A,[vl>(y),vl/'(x)] =6ix-y)=5{y-x) 



(5.18) 



which is the distributional differential equation for the fundamental solution of the 
Laplacian. This is straightforward to solve and yields 



[^ix),^'iy)] 



rO=yO /^Tiz"^ \x — y\ 

such that the result of lemma 13.11 can be readily applied 

[x{x),x{y)] „ = [x'ix),x'{y)] 

X — y 







i iey} 1 

[x{x),x'{y)] o = —^(^~y)~-A — 2^ — ^ 

x"=y" Ana^ \x — y\ 



This proves the claim. 



(5.19) 

(5.20) 
□ 
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This result has intriguing imphcations, as \1/ and x have a clear physical inter- 
pretation but the "natural" quantization procedure suggests that they are non-local 
fields. Such a situation may point to an underlying non-commutativity of space- 
time, as the cause of this non-locality. It is worth noting that "non- commutative 
infiation", i.e. infiationary cosmological scenarios based on a non-commutativity of 
spacetime, is a topic of evolving interest during the past years, see e.g. PQ [S] [12] [21] • 
The above result may be interpreted as indicating that it is "natural" to assume 
spacetime non-commutativity in infiationary scenarios. 

6 Weyl Quantization 

In this section the algebraic Weyl quantization procedure is applied to the dynamic 
system of perturbations. As a general reference for algebraic methods in quantum 
field theory see [T7j, for the application of Weyl quantization in quantum field theory 
on curved spacetime see [21] • An application of the Weyl algebra in the context of 
loop quantum gravity is presented in [2S]- The Weyl quantization procedure has the 
benefit that the operators involved are bounded if represented on a Hilbert space, 
which makes it possible to define a C* norm on the Weyl algebra. 

For the quantization procedure which is applied here it is of great importance that 
the constraints are second class constraints. For systems with first class constraints 
quantization is much more complicated [16] and methods such as refined algebraic 
quantization [2] and more recently a functional approach to the Batalin-Vilkovisky 
formalism [2] [12] have been developed. For the case at hand a rather conventional 
approach suffices. 

The general pattern that will be followed is the same as in the previous subsec- 
tions. First the system without constraint!! is considered and then the impact of 
the constraint on this system is investigated. The present subsection does not pro- 
vide additional information on the dynamical system but clarifies the impact of the 
constraint at the level of algebraic quantum field theory. 

As we will deal with the quantum system here and not the classical system, the 
question whether quantization commutes with reduction comes into play. In the 
present setting, reduction on the classical level can be split up into two parts. First 
the Poisson bracket on phase space is replaced by the Dirac bracket, which induces a 
presymplectic form, and in the second step a projection to the symplectic subspace 
is performed. It is not clear how the first part of the reduction procedure should be 
implemented at the quantum level, while the second part of the reduction procedure 
gives rise to a well defined algebraic procedure, such that quantization commutes 
with the second part of reduction by construction. 

Starting out from the dynamic system for \E' and x is unsuitable for a consideration 
from the point of view of algebraic quantum field theory as the system contains cou- 

^We speak only of one constraint here, as the second constraint is derived from this constraint 
as a secondary constraint. 
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pling of the fields, which undermines a completely rigorous treatment with current 
methods. Instead it would be preferable to start from a system of two independent 
free scalar fields. To achieve this, it is necessary to apply the constraints to arrive 

2 

at a dynamic system for ^' = :^\E' and u. 

*"-A*+ (^-2e + 6 + e6 -6^ - ^H'^^ = (6.1) 

u"- Au+ (^-2 -2e + 3S- 2e^ + 3e5 - + ^ j "^^^ = (6-2) 
and the constraint takes the form 

^' + {1 + e - 6)71^ = —u = 4nGu (6.3) 



Ignoring the constraint, ^ and u are free scalar fields and the symplectic space 
3^ of their modes is straightforward. Now in the present treatment the reduction 
of phase space is done in two steps. First the Poisson bracket is replaced by the 
degenerate Dirac bracket, leading to a space ^ with a symplectic subspace. Then 
a projection to the maximal symplectic subspace is performed by factoring out the 
primary and secondary constraint, leading to a symplectic space ^ jC 

To tackle the problem of quantization, it is useful to take the spatial Fourier 
transform of the fields, as this allows to make direct contact with the previous 
treatment. In this context, the Dirac form is defined as in equation 04.81) . where 
the mode Dirac form is of course different from the one given in equation 04.71) . as 
we are working with different coordinates on the symplectic space here. However, 
it can be expected that the limit A; — ?■ should be of the same order of divergence 
as in equation 04.71) because the coordinates are related by a /c-independent linear 
transformation. Therefore the integral can be expected to exist. The constraint in 
mode form is 

Ci,fc := *1 + (1 + £ - (5)OTfe - 47rGMfc = 
and the Dirac bracket of the constraint with an arbitrary initial value set is 



i ^ \ 
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The algebras of observables for ^ and u, the Weyl algebras Ql-^f/u, can be con- 
structed in the standard way. The algebra of the full system is then given as the 
tensor product algebra 21 = 21^ ® 21^. Now the constraint cannot be simply inter- 
preted as a relation in this algebra that has to be factored out to restrict to the 
constrained system, because it is incompatible with the commutation relations. As 
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discussed above, applying the constraint is not very complicated at the level of field 
modes; at the level of the Weyl algebra the Weyl relations have to be replaced with 
the modified relations derived from the mode Dirac forms, profoundly changing the 
algebra. 

For the thus acquired algebra Qi^ corresponding to ^ the constraint defines a 
non-trivial ideal, as the Dirac form is constructed to be degenerate in exactly this 
sense. This is in contrast to 21 which, as a Weyl-Algebra corresponding to a non- 
degenerate Poisson form, is simple. The Weyl operators are defined with respect to 
a set of mode initial values 

and the Dirac- Weyl relation is 

W{fi, /2, /a, U)W{gi, c/2, 93, Qi) 



exp 



2] h ' 93 

\ [\Ul \94 / }J 



W{h + 91, f2 + 92, h + 93, h + 94) (6.5) 



This makes it possible to formulate the constraint at the level of the Dirac- Weyl 
algebra 21^ as 



V/ e ^(I 



W{f,{l + e-S)'Hf,0,-4nGf) = l 



where ^(M+) denotes the space of Schwartz functions. The constraint is compat- 
ible with the Dirac- Weyl relation 06. 5p by virtue of equation 06.41) . The secondary 
constraint can be calculated and cast into a similar form at the level of the alge- 
bra. The constraint relations define two ideals on the algebra and are factored out. 
The resultant algebra is denoted as 2ts?/C' and corresponds by construction to the 
reduced phase space ^/C. This is illustrated in the following diagram, where it 
should be noted that the connection between 21 and 21^ is not clear at the purely 
algebraic level. Instead 21^ is constructed using results at the classical level. 



Dirac 



Reduction 



21 



Quantisation 

7 



-> 21 



Quantisation 

Reduction 
& 



^/C 

Quantisation 



The commutation of quantization and reduction has been investigated in a plethora 
of models with a variety of tools. For a precedent to the type of commutative di- 
agram presented here see e.g. j23] where a gauge theoretic model with first class 
constraints is investigated applying much more sophisticated technical tools than 
are used in the present work. 
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States from the state space S of the tensor product algebra 21 can be freely chosen, 
for instance it is possible to pick states which have the Hadamard property for both 
fields. The simplest states of this kind are simply tensor products of Hadamard 
states for the single field algebras 2l*/„. When the constraint C is applied, the 
resultant algebra 2t^/C' has a very different state space Sc- In the construction 
applied here, the isomorphism 2t^/C' ~ 2l„ is implied, which permits a canonical 
injective, unit preserving C*-homomorphism a : Sl^^/C — 21. The dual map to this 
homomorphism is a positive map a' : S Sc which means that states on 21 give 
rise to states on 2ls?/C' by simply restricting to their action to 21^ C 21. 

The explicit form of the relations between states in S and states in Sc is greatly 
simplified by our choice of dynamical variables. If one chooses different dynamical 
variables, like \1/ and the construction of a corresponding algebra 21' is much more 
troublesome and a simple relation between 21' and 2t'^/C' ~ 21^/C' need not exist. 

7 Conclusion 

We have shown that the fact that the dynamic system of linear scalar cosmological 
perturbations has only one degree of freedom makes it impossible to find a quan- 
tization procedure such that the metric and field perturbation both satisfy canoni- 
cal commutation relations. By means of the Dirac formalism for the treatment of 
dynamic systems with constraints the standard literature choice for the canonical 
variables was justified from purely formal arguments. The Weyl quantization of the 
system is tractable and leads to a generic scalar field algebra. The physical degrees 
of freedom, namely the metric and field perturbations emerge as non-local fields 
from the quantization procedure. 
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